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ABSTRACT 

The possible formation of Direct Collapse Black Holes (DCBHs) in the first metal-free atomic 
cooling halos at high redshifts (z > 10) is nowadays object of intense study and several meth¬ 
ods to prove their existence are currently under development. The abrupt collapse of a massive 
(~ 10^ — 10^ Mq) and rotating object is a powerful source of gravitational waves emission. 
In this work, we employ modem waveforms and the improved knowledge on the DCBHs 
formation rate to estimate the gravitational signal emitted by these sources at cosmological 
distances. Their formation rate is very high 10^ yr“^ up to z ~ 20), but due to a short du¬ 
ration of the collapse event (~ 2 — 30 s, depending on the DCBH mass) the integrated signal 
from these sources is characterized by a very low duty-cycle (V ^ 10“^), i.e. a shot-noise 
signal. Our results show that the estimated signal lies above the foreseen sensitivity of the 
Ultimate-DECIGO observatory in the frequency range (0.8 — 300) mHz, with a peak ampli¬ 
tude Hgu, = 1.1 X 10“^'* at Vmax = 0.9 mHz and a peak Signal-to-Noise Ratio SNR ^ 22 at 
1 / = 20 mHz. This amplitude is lower than the Galactic confusion noise, generated by binary 
systems of compact objects in the same frequency band. For this reason, advanced techniques 
will be required to separate this signal from background and foreground noise components. As 
a proof-of-concept, we conclude by proposing a simple method, based on the auto-correlation 
function, to recognize the presence of a D ^ 1 signal buried into the continuous noise. 

The aim of this work is to test the existence of a large population of high-z DCBHs, by ob¬ 
serving the gravitational waves emitted during their infancy. 

Key words: gravitation - waves - cosmology; theory, early Universe - galaxies; high redshift, 
formation - quasars; supermassive black holes 


1 INTRODUCTION 

Gravitational Wave (GW) observations are likely going to open a 
new window on the Universe within this decade. The different emit¬ 
ters of electromagnetic radiation and GWs are responsible for the 
dissimilar spatial regions they convey information about. The elec¬ 
tromagnetic radiation is emitted by charged particles with relatively 
small wavelengths, carrying information about the physical con¬ 
ditions of small spatial regions. On the contrary, the entire mass 
and momentum distributions of an astronomical source produce the 
emission of GWs, which carry information about large spatial re¬ 
gions, with very long wavelengths. Moreover, electromagnetic ra¬ 
diation and GWs couple very differently with the matter; while the 
former is easily absorbed or scattered, the coupling of the latter is 
extremely weak. This effect allows the GWs to travel cosmological 
distances without being affected, but also makes them very hard to 
be detected (see e.g. |Sathyaprakash & Schutz|2009^ . 

First-generation GWs detectors (e.g. LIGO, Virgo) have been 
successfully commissioned, and the development of the next ad¬ 
vanced sensitivity ground-based detectors (Advanced LIGO, Ad¬ 
vanced Virgo) is well underway. However, the construction of the 


planned space-based detectors (DECIGO, eLISA) is still in the de¬ 
velopment phase (see e.g. |Kawamura et al.|2011 [ |eLISA Consor-| 
|tium et al.|2013k 

The GW spectrum is conventionally divided into several re¬ 
gions, from extremely low (~ 10“^® Hz) to high (~ Hz) fre¬ 
quencies. In the low frequency range (approximately from 10“^ Hz 
to 1 Hz, see |Camp & Comish||2004| for an extensive review) the 
most important sources are Galactic compact binaries, massive 
Black Holes (BHs) binary mergers, massive BHs capture of com¬ 
pact objects, the collapse of super-massive stars and the primordial 
GWs background (see |Amaro-Seoane et al.|2013| for an extensive 
review on the subject). 


In particular, the collapse of Super-Massive Stars (SMSs), the 
first massive objects formed in unpolluted atomic cooling halos at 
z > 10 may lead to a significant emission of GWs, depending on 
the collapse asymmetry. A possible mechanism for a highly asym¬ 
metrical collapse is the development of a dynamical bar-mode in¬ 
stability as the super-massive star cools (Smith et^r]|1996[ |Saijo| 


|et al.|2002[|Shapiro|2003[[Franci et al.|2oi3[ ). This may be likely if 

viscosity and magnetic fields are insufficient to keep the star rotat 


© 0000 RAS 
















2 F. Pacucci et al. 


ing uniformly during the cooling phase. Given enough energy and 
a long enough lifetime of the bar mode, a significant fraction of 
the rest energy of the star could be lost as gravitational radiation 
(see e.g. [Schneider et al.|2000| l. If this signal is actually observable 
by current or future GWs surveys, it would provide a highly valu¬ 
able tool to prove the existence of Intermediate Mass Black Holes 
(IMBHs) at high redshifts. 

In the commonly accepted ACDM cosmological scenario, 
the formation of the first stars and black holes occurred at z = 
20 — 30 ( |Miralda-Es c ude|2003[|Bromm & Yoshida|20ir||Volon-| 
|teri & Bellovary|2012a[ ) in molecular cooling halos, i.e. dark mat¬ 
ter structures with virial temperatures below 10"^ K and masses be¬ 
low Mh ~ 10® Mq. These structures cooled their metal-free gas 
through molecular hydrogen line emission. Under intense UV ir¬ 
radiation, particularly in the Lyman-Werner (LW) band (11.2 — 
13.6 eV), the molecular hydrogen is photo-dissociated, so that the 


cooling is quenched (see [Machacek et al.||2001[ [Fialkov et al.| 
|2012^ . On the contrary, larger metal-free halos (with masses above 
~ 10® Mq and temperatures T^ir > 10"^ K), when illuminated 


by LW photons with fluxes above a threshold J* dLoeb & Rasio 

1994| Eisenstein & Loeb|[1995| Begelman et al.|[2006| [Lodato & 

Natarajan|2006 

[[Regan & Haehnelt|2009[ Shang et al.|2010 John- 

son et al.|2012 

[Agarwal et al.|2012[[Latif et al.|2013[[Inayoshi & 

Omukai|2012| 

Sugimura et al.|2014[ Dijkstra et al.|2014l are the 


ideal environment for the formation of Direct Collapse Black Holes 


(DCBHs). Several theoretical works l|Bromm & Loeb|2003 

Begel- 

man et al.|2006 Volonteri et al.|2008[ Shang et al.|2010[[. 

ohnson 

et al.|2012[ [Ferrara et al.|2014|l have shown that the result of this 


collapse is the formation of a DCBH of mass 10^ — 10® Mq . The 
collapse of SMSs is mainly driven by General Relativity (GR) insta¬ 


bilities, as described in |Shibata & Shapiro| ( |2002l l. The subsequent 
accretion of mass contributes to the final mass of the IMBHs. 


The existence of IMBHs at high-redshifts would provide a 
very convenient solution to two problems still unsolved by mod¬ 
ern Cosmology, namely: (i) What are the sources responsible for 
the Cosmic Infrared Background jCappelluti et al.||2013| |Yue| 
|et al.|20I3| l? (ii) How did the current Super-Massive Black Holes 
(SMBHs) form (see |Volonteri & Bellovary|2012b| for a recent re¬ 
view and references therein)? We propose here a tentative test of 
the IMBH hypothesis based on the detection of the gravitational 
waves (GWs) emitted during the collapse phase. 


Several works have been focused on the theoretical determi¬ 
nation of the GWs waveform for the most important sources in 


the low frequency range, starting from the seminal paper Saenz 
& Shapiro ([jjTsjl. For example, Schneider^etayj2000 Ott et al. 


[Sekiguchi & Shibata||2005[ |Ott| 2009[ Li & Benacquista 
have developed waveforms for the asymmetrical collapse of a 
SMS into a black hole both with simple theoretical arguments and, 
more recently, with numerical simulations. For instance, [Reisswi ^ 
|et al.|p0T3t have investigated a non-axisymmetric instability which 
leads to the formation of a binary black hole system within the col¬ 
lapsing SMS. Employing this modem waveform for the collapse of 
DCBHs at high redshifts (z ~ 15) we study their observability with 
the aid of recent estimates of their formation rates. 


2004 


2010 


The outline of this paper is as follows. In §2 we present 
an overview of the theory describing the GWs emission, spec¬ 
ifying the waveform used. In §3 we describe the results and 
the possible observability of these gravitational signals. Finally, 
in §4 we state our conclusions. Throughout, we adopt recent 
Planck cosmological parameters ( [Planck Collaboration et al.|2013^ : 
{Urn, Ua, Ub, h, rig, (Ts) = (0.32,0.68,0.05, 0.67,0.96,0.83), 


2 THEORETICAL BACKGROUND 

The output of a GWs detector is a time series s{t) which includes 
the instrament noise n{t) and the response to the gravitational sig¬ 
nal h{t): 


s{t) = P+{t)h+{t) + Px{t)hx{t) + n{t) (1) 

The instrument response is a convolution of the antenna patterns 
P+{t) and Px{i) with the two GW polarizations h+ and hx- The 
antenna patterns depend on the sky location and on the emission 
frequency of the source and they are simple quadmpoles in the case 
of emission wavelengths which are large when compared to the de¬ 
tector baseline. The signal analysis is usually performed in the fre¬ 
quency domain since, in this representation, the noise is usually 
assumed to be uncorrelated and gaussian in each frequency bin. 

The information contained in the time series is generally rep¬ 
resented in the Fourier domain as a strain amplitude spectral den¬ 
sity, h{v). This quantity is defined in terms of the power spectral 
density: 


Ssiiy) = s*{v)s{v) (2) 

The tilde operator is the Fourier transform of the time series: 

/ + 00 

df. (3) 

-oo 

and the star indicates the conjugation operation in complex num¬ 
bers. The power spectral density Ss {v) has units of time or, as usu¬ 
ally indicated, Hz“^. The strain amplitude spectral density is then 
defined as: 


h{v) = a/ Ssiu) (4) 

This quantity has units and is the one commonly reported 

in plots showing gravitational signals. From the prediction of the 
time evolution of the gravitational signal shape h{t), derived from 
the instrumental time series s{t), it is then possible to compute the 
strain amplitude spectral density h{v). 

It is standard practice to quote the strength of a gravitational 
signal in terms of the energy density per logarithmic frequency in¬ 
terval, dpgm /d\nv, scaled by the energy density needed to close the 
Universe: 


Ugro{l^'} — 

Peril 


dpgw 

dlni' 




V^Ss{v) 


(5) 


where Hq is the Hubble constant, perit is the critical density and 
pgru is the energy density in a GW, given by: 


Pg-w — 


1 

32)(G 


<h > 


( 6 ) 


where G is the gravitational constant and the brackets denote a spa¬ 
tial average over the wavelengths. 

The actual prediction of the shape h{t) is far from being trivial 
and strongly depends on the nature of the source and on its physi¬ 
cal properties. Only recently a number of accurate waveforms have 
been proposed for binary systems and for collapsing SMSs (see the 


works Schneider et aL|2000[ 

Li & Benacquista|2010[|Li et al.|2012| 

[AJith et al.[[2011 [Pan et al. 

2014 1 . The form of the gravitational 


waveform template plays a crucial role for data analysis: very often, 
only a detailed prediction of the signal shape allows to disentangle 
it from the receiver noise and/or from other sources of gravitational 
radiation. We will discuss this point more carefully in the following 
subsections. 
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Gravitational Waves from Direct Collapse Black Holes Formation 3 


2.1 Waveform for IMBH Collapse and Ringdown 


A DCBH can emit GWs at its formation through the asymmetric 
core collapse of its progenitor. Indeed, in order to emit gravitational 
radiation, a physical system needs to be non-spherical: one possibil¬ 
ity for having a non-spherical collapse is requesting the primordial 
halo, hosting the to-be-formed BH, to be in rotation. DCBHs are 
formed during a brief era of cosmic time (13 < « < 20, Yue et al. 
| 2014^ due to the collapse of the inner part of metal-free halos. Ihe 
angular momentum of high-redshift halos is small ( |Davis & Natara-| 
[jan|2010) and is preserved during the collapse. The newly formed 
DCBH is then characterized by low values of angular momentum, 
which is later increased by accretion and merging events. 

We focus on BH progenitors, i.e. the inner section of the col¬ 
lapsing halo, with mass M, mean mass density p and dimensionless 
spin parameter a defined as: 


Jc 

GAF 


(7) 


where J is the angular momentum and c is the speed of light. Note 
that a — Q denotes a non-rotating BH, while a = 1 is a maximally 
rotating Kerr BH, although we note that |Thom^ ( |1974^ showed that 
accretion driven spin is limited to a = 0.998. Magnetic fields con¬ 
necting material in the disk and the plunging region may further re¬ 
duce the equilibrium spin: magneto-hydrodynamic simulations for 
a series of thick accretion disks suggest an asymptotic equilibrium 
spin at a ~ 0.9 ( |Gammie et al.||^04| ). Rotating BH progenitors 
will distort to an oblate spheroid shape despite their own immense 
gravity. The collapse, under this non-spherical geometry, causes the 
emission of a huge amount of energy through GWs. The final evo¬ 
lution of the collapsing halo is usually characterized by two distinct 
phases: 


(i) Collapse: under the influence of its overwhelming self¬ 
gravity, the BH progenitor collapses in a time scale equal to the 
dynamical time tdyn ~ 1 / fGp. 

(ii) Ringdown: the collapsed non-spherical object undergoes 
strong oscillations which progressively allow it to acquire a spheri¬ 
cal shape. 


As we shall see, the maximum of gravitational radiation is released 
at the end of the collapse, when the BH bounces back and starts the 
subsequent ringdown phase. 

In the work by |Li & Benacquist^ j2010^ , the core collapse 
of the SMS is approximated as an axisymmetric Newtonian free- 
fall of a rotating relativistic degenerate iron core. In addition, the 
collapse waveform is reasonably well modeled by an exponential 
growth. Following this paper, the time-series of the waveform for 
the ringdown phase is expressed analytically as: 


hRoit) = A——eyip{—-RVot/Q)cos{2'KVot) (8) 
c-^dL 

where dL is the luminosity distance between the source and the 
detector, Q = 2(1 — g{a) = 1 — 0.63(1 — vq = 

Fg{a)/{2 tvGM). Note that is the typical timescale of GWs 

emission. In addition: 

A = y| (i + »(«)-« (9) 

where e is the fraction of the initial mass of the BH which is trans¬ 
formed into GWs radiation. Numerical simulations have shown that 
approximately 1% of the final BH mass is emitted in GWs so, 
throughout, we adopt the value e = 0.01 = 1% (see | |Buonanno| 
let al.|200'7l [Abbott et al.|200^ for details). In the latter work, the 


authors investigate the presence of possible GWs burst signals in 
the high-frequency range 1 — 6 kHz, without finding any evidence 
of them, but putting an indirect upper limit on the emitted energy. 
Although the survey was dedicated to a different frequency range, 
we make use of the same upper limit, in accordance with |Li & Be-| 
|nacquista|j 2010 ^ . 

The waveform for the collapse is found by analytically fitting 
a numerical relativistic simulation in |Li & Benacquist^j2010^ . The 
time series has the following general form: 

hc{t) ~ a + (10) 

P 

where a is the spin parameter, M is the mass, p is the mean mass 
density of the collapsing object, and a and 7 are two free param¬ 
eters that we determine by imposing that hc(t = 0) — 0 and 
hc{t = tc) = hnoit = tc), i.e. the gravitational signal is null at 
the initial time and the two waveforms match at the collapse time 
tc. The collapse time is found by imposing its identity with the dy¬ 
namical time of a quasi-spherical object with a mean mass density 
P- 


tc 


tdy 


1 

y/Gp 


( 11 ) 


Calling V the value of the ringdown waveform at the collapse time 
tc, our final form for the collapse waveform is the following: 


with 


he it) = ^ - 1) 

7= 

Mtc \a 1 


( 12 ) 


(13) 


The physical parameters in the model are: the mean mass density 
of the collapsing object p, the BH mass M and the spin parameter 
a. The range of their values used in our calculations is detailed in 
Sec.[^ Fig.[^shows an example of time series for the gravitational 
strain: the collapse and subsequent ringdown phases are evident and 
separated by the vertical line at tc ~ tdyn ~ 0.7 s. Fig. shows 
the spectral strain for a single source with the parameters reported 
in the caption. The peak amplitude ~ Hz“^/^ is reached at 

i/ ~ 3 X 10“^ Hz. 

Other waveforms have been proposed for the collapse of SMSs 
with lower masses. For example, |Suwa et^ ( |2007) > employed GR 
numerical simulations and the standard quadrupole formula de¬ 
veloped in [Moenchmeyer et al.| jl991^ to study the collapse of 
~ 500 Mq SMSs, obtaining a peak power in the range 10 — 100 Hz. 
Due to important differences in the collapse modeling, a direct 
comparison between the classical quadrupole formula derived in 
[Moenchmeyer et ^ ( [1991^ and the one used in the present work is 
hard to obtain. 


2.2 Other Sources of Gravitational Signals in the Low 
Frequency Band 

Most of the gravitational power emitted from the DCBH collapse 
falls in the low frequency range, approximately from lO”'^ Hz to 
IHz (see [Camp & Cornish|2004| for an extensive review). In this 
range the most important sources of GWs are Galactic and ex- 
tragalactic compact binaries, massive BHs binary mergers, mas¬ 
sive BHs capture of compact objects and the primordial GW back¬ 
ground. The last one is too weak to be of any concern in this work 
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Figure 1. Gravitational strain as a function of time for the collapse of a 
10® Mq BH at 2 = 15, with a mean mass density p = 10^ g cm“®. The 
separation between the collapse and the ringdown phases is mai'ked by the 
veifical blue line. 



Figure 2. Spectral strain for a single source, due to the collapse of a 10® Mq 
BH at 2 = 14, with a mean mass density of p = 10^ g cm“®. The power 
peak is reached at t/ ~ 3 X 10“® Hz. The sensitivity curve for the DECIGO 
obseiwatory is reported as a shaded orange line. 


(see e.g. |Yagi & Seto|2011| >. This is not the case with the Galactic 
and extragalactic compact binaries. 

The unresolved background produced by the GWs emission 
of Galactic and extragalactic compact binaries acts as a confusion 
noise in this frequency rang 


ney|20031 INelemans 

20091 ISathyaprakash & Schutz|2009||Regim- 

bau & Hughes|2010 

IMarassi et al.|2011 |Regimbau et aL|2012J. 


This is produced by the coalescence of two neutron stars (NS-NS), 
two black holes (BH-BH) or a neutron star and a black hole (NS- 
BH) in our Galaxy or in the 2 < 5 Universe. The duty-cycle (i.e. 
the percentage of one period in which a signal is active, see a more 
complete definition below) for these events, for future space-bome 
gravitational observatories like DECIGO, will be higher than unity, 
i.e. the signals overlap and create a confusion noise of unresolved 
sources. 

The merging event for SMBHs also produces very energetic 
GWs ( [Kocsis & Loeb|2008| > that may represent an important fore¬ 
ground contribution to the detected gravitational signal. One very 
simple method to reproduce the gravitational signals generated by 


SMBHs merging events is presented in |Li et aT] ( |2012^ where the 
interested reader is referenced to for a detailed description. How¬ 
ever, the signal generated by these events is very different (in terms 
of both the duty-cycle and the characteristic energy) from the one 
of our interest and their separation should be simple. 


2.3 Formation Rates of DCBHs 

We employ the mass density of high-redshift DCBHs predicted in 
the work |Yue et al.| ( [^13| >, where they present the rate of formation 
per unit volume of DCBHs, ip{z), that accounts for the observed 
Cosmic Infrared Background (CIB) fluctuation excess. The Uni¬ 
verse enters the DCBHs era at 2 ~ 20 when a large fraction of 
atomic cooling haloes are experiencing DCBHs formation. Their 
formation is suppressed after 2 « 13, so that the DCBHs era lasts 
only ~ 150 Myr of the cosmic history (see |Yue et al.|201^ for de¬ 
tails). Then, the number of GWs sources formed per unit time out 
to a given redshift 2 can be computed by integrating the cosmic 
DCBHs formation rate per unit volume <^( 2 ), see Fig. ^ 

It is important to note that the |Yue et'aL]j2013||2014| l estimated 
DCBH high -2 mass density (~ 2 x 10® Mq Mpc“®) exceeds the 

|Tremaine||2002| computed through the standard Soltan argument 
(which derives this quantity from the quasar output, see |Soltan|1982| 
for the original paper). This point is worth some discussion. 

First, the |Yue et al.| ( |2013||2014) DCBH abundance is required 
in the extreme case in which most of the CIB fluctuations in the 
range (3.6 — 4.5) /rm are to be explained by DCBHs. This might 
not be the case, and the difference between the local BH density 
and the high -2 DCBH one might be eased. 

Second, even if the CIB constraint is imposed, the BH mass 
density adopted here cannot be ruled out by current observations. In 
fact, estimates derived from SMBH demographics refer to BHs that 
(i) are relatively massive, and (ii) reside in structures that are iden¬ 
tified as “galaxies”. As for point (i), there could be a large number 
of BHs with masses ~ 10^“® Mq that are too faint to contribute 
sensibly to the X-Ray Background (XRB) and/or are very ineffi¬ 
ciently accreting (see e.g. studies considering slim-disk solutions 
like |Begelman & Meier|[T98^ |Paczynski & Abramowicz||1982[ 
Sadowski||2009| |2011| and the recent work by |Volonteri & Silk| 
20 14^ . These objects could be easily misidentified with other X- 
ray sources, e.g. High-Mass X-Ray Binaries (HMXBs) and Ultra- 
Luminous X-ray sources (ULXs). The additional possibility (ii) ex¬ 
ists that a considerable fraction of DCBHs resides in baryon-free 
dark matter halos, as suggested by recent studies, e.g. |Pacucci &| 
|FetTara|j2015| l. These authors have simulated the accretion process 
onto a DCBH of initial mass 10® Mq at 2 ~ 10, embedded in a 
dark matter halo with a gas content of ~ 10^ Mq, finding that after 
~ 142 Myr about 90% of the initial gas mass of the halo is accreted 
onto the compact object. If these objects are not embedded in larger 
gaseous structures they could completely escape our census. 

Finally, [Kormendy & Ho| ( |2013| ) suggested a revision of the 
BH mass scaling relations, indicating that the local mass density 
in black holes should be increased by up to a factor of five with 
respect to previously determined values. As already pointed out, a 
possibility to explain a large BH mass density is that most of their 
growth occurs via radiatively-inefficient channels. In a recent paper, 
IComastri et al.|(20I5[ l showed that it is possible to accommodate a 
large fraction of heavily buried, Compton-thick AGNs, without vi¬ 
olating the limit imposed by the hard X-ray and mid-infrared back¬ 
grounds spectral energy density. 

Thus, before we can firmly exclude the existence of a sub- 


present-day SMBH mass density (~ 2 x 10® Mq Mpc 


Yu & 
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Figure 3. The event rate and the duty-cycle for the formation of DCBHs as a 
function of redshift. The upper lines (solid and dotted ones) are cumulative, 
i.e. obtained summing all the values up to the given redshift. 


stantial population of intermediate mass black holes at high -2 one 
would need to predict in detail where these DCBHs end up locally, 
what is their accretion history, and what are the best strategies to 
detect them. The aim of this work is to explore the detectability 
of these objects via their gravitational wave signal, at least in the 
optimistic, but physically motivated, case considered here. 


2.4 Calculation of Event Rates, Duty-Cycles and S/N Ratios 


An important parameter to describe the signal in the time domain 
is the duty-cycle, T). This is defined as the ratio between the typi¬ 
cal duration of the signal emitted by a single source and the aver¬ 
age time interval between two successive emissions. When D > 1, 
the overall signal is continuous (i.e. the different signals overlap in 
time), conversely if D < 1 the resulting background is character¬ 
ized by a shot-noise structure. From the value of iplz) and from the 
Initial Mass Function (IMF) <E>(M) of the population, it is possible 
to compute the duty-cycle as: 


dVjz) 

dz 


(p{z) dV (1 + z) 
(1 + z) dz i/Q 



(14) 


Here , the IMF of IMBHs seeds (normalized by the mass) is taken 
from 


Ferrara et al. 

emission (see Sec. The double appearance of the term (1 + 
to take into account the redshift dependence of both vq and (p. 


1 2014 I and Vq is the typical timescale of GWs 

z) is 


The Signal-to-Noise Ratio allows a direct comparison between 
the GW amplitude and the instrumental sensitivity, assessing what 
sort of sources will be observable against noise. A stationary noise 
is usually described by its Power Spectral Density (PSD): as a con¬ 
sequence, these comparisons are most commonly made in the fre¬ 
quency domain. 

In this work, we deal with short-lived signals, which have wide 
bandwidths and the actual observation time of the source is not rel¬ 
evant in the determination of the Signal-to-Noise Ratio (SNR). It 
is useful to define the dimensionless noise power per logarithmic 
bandwidth: 

hl{p) = (15) 

where hn{v) is called the effective GW noise. From the signal side, 
we may define the characteristic signal amplitude: 

hc = v\h{v)\ (16) 

which is dimensionless. This last quantity is to be compared with 
the effective GW noise h„ (v) to obtain a rough estimate of the SNR 
of the signal: 

SNR{v) = (17) 

hn{v) 

The sensitivity curve for Ultimate-DECIGO has been taken from 
|Marassi et al.| ( |200^ and |Yagi & Seto|j2011| >. 

The SNR with respect to the instrumental noise is not the only 
parameter to be taken into account in order to assess the observ¬ 
ability of a source. The signal also needs to be distinguished from 
all the other components received in the same frequency band, as 
described previously in this Section. The observability will be ad¬ 
dressed more carefully in Sec.[^ 


3 RESULTS 

In this Section we present the results of our calculations and we 
discuss on the observability of the signal. In Sec.j^we presented the 
employed waveform for the gravitational signal generated by the 
collapse of a high-redshift atomic cooling halo, with the subsequent 
formation of a DCBH. The waveform is a function of the mean mass 
density p, the total mass M and the spin parameter a (see Eq.|12[>. 
The amplitude of the detected gravitational signal depends also on 
the luminosity distance from the source, i.e. from its redshift. 

In order to simulate the gravitational signal generated by a re¬ 
alistic ensemble of the population of high-z DCBHs, we let these 
parameters vary into different ranges of values, summarized in the 
following Table 1. 


The event rate and the duty-cycle are shown in Fig. The 
event rate for all-sky observations is very high, reaching a peak of 
~ at 2 ~ 22. However, due to a rather short duration of 

the GW burst (~ 2 — 30 s depending on the DCBH mass, see Fig.[TJ 
the duty-cycle is very small, of order 10~^. This result suggests that 
the gravitational signal produced by the collapse of high-redshift 
DCBHs is a shot-noise. The decrease, with increasing redshift, of 
both the event rate and the duty-cycle (and the corresponding level¬ 
ing of their cumulative values) is ascribed to the decrease in p{z), 
as in|Yue et al.|(|2013|(. Both the event rate and the duty-cycle are 
comparable with the results in |Marassi et al.|j2009j who, however, 
considered the contribution to the stochastic background of Pop III 
stars with masses in the range (300 — 500) Mq that collapse to 
BHs. 



MIN 

MAX 

Bins 

z 

14 

22 

10 

Logp[g cm"®] 

6 

8 

10 

Log M [Mq] 

4.5 

5.7 

10 


The redshift range has been chosen in accordance with the results 
of|Yue et^j2013t, while the mass range to be compatible with the 


IMF values in Ferrara e t al.H201 4^. The den sity interval has been 
taken from simulations in Li & Benacquista| ( |2010|l, while th e vari¬ 
ation of the spin parameter a for high-z BHs from |Volonteri|p(jT^ , 
|Volonteri et al.| ( [?013| >. Generally speaking, the rotational energy of 
primordial BHs is very low compared with local ones: the spin is 
progressively increased by the accretion of matter onto the compact 
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Figure 4. Effect on the strain amplitude h(u) of the BH mass and of its red- 
shift. In this case the calculation has been performed in narrow ranges of the 
parameters, as indicated in the legend. The orange region is the sensitivity 
for the future Ultimate-DECIGO observatory. 



Figure 5. As in Fig.[^ but for the spin parameter and the mean density. In 
this case, the single values used for the parameters are shown in the legend. 

object and by merging events. For this reason, our average value for 
the spin parameter is a ~ 0.05 (see also |Davis & Natarajan|201^ . 

The variation of these physical quantities have different im¬ 
pacts on the produced gravitational signal. A study of their effect 
has been sketched in the following Fig. I^and Fig. 0 In these fig¬ 
ures the range of variation for the different parameters is restricted 
to the values shown in the legend. 


The effect of the redshift on the gravitational signal is very 
small. This is due to the fact that the whole range of redshift val¬ 
ues spans only ~ 150 Myr of cosmic history. For higher redshifts, 
the observed gravitational signal is shifted to lower frequencies, as 
expected. On the contrary, the effect of the mass is greatly impor¬ 
tant: lower masses correspond to higher frequencies of the peak 
{vo oc 1/M, see Sec.[^, but also to an overall amplitude of the sig¬ 
nals ~ 2 orders of magnitude lower, at the extreme detection limit 
with the DECIGO observatory. The spin parameter affects mainly 
the shape of the spectral strain amplitude and its peak frequency. 



Figure 6. Gravitational strain h(v) for the total ensemble. The orange re¬ 
gion is the sensitivity for the future Ultimate-DECIGO observatory. The 
dashed lines indicate the Galactic Confusion Background (GCB) due to the 
GWs emission from double neutron stars (NS/NS), which is the dominant 
foreground component (see |Yagi & Seto|2011^ . 

which is shifted to higher values for larger values of a, while the 
amplitude is not affected. The mean mass density p has virtually 
no effect on the peak frequency and on its amplitude, but varies 
the shape of the spectrum at higher frequencies. This is due to the 
fact that the mean mass density directly controls the duration of the 
collapse phase, tc ~ tdyn ~ 1/^Gp, i.e. the transition from the 
collapse to the ringdown phases. 

The gravitational strain h{i') for the total ensemble is reported 
in Fig. 1^ The orange region is the sensitivity curve of Ultimate- 
DECIGO, as obtained from [Marassi et al.| and |Yagi & | 

|Seto| pOl 1| >, and it has been used to compute the SNR, reported 
as a solid red line in Fig. The estimated signal lies above the 
foreseen sensitivity of Ultimate-DECIGO in the frequency range 
(0.8 — 300) mHz, with a peak amplitude 0^™ = 1.1 x 10“®“* 
at Vmax = 0.9 mHz. The Signal-to-Noise ratio reaches a maxi¬ 
mum value of SNR ~ 22 at= 20 mHz. A source is considered 
detectable if the resulting SNR exceeds some standard threshold 
value, typically between 5 and 10 ( [Plowman et al.|20T0| , so the sig¬ 
nal generated by the collapse of high-redshift atomic cooling halos 
into DCBHs is detectable by the future Ultimate-DECIGO obser¬ 
vatory. 


However, in this frequency interval the sensitivity of Ultimate- 
DECIGO is limited by the unresolved background produced by the 
GW emission of Galactic compact binaries that acts as a confu¬ 
sion noise (|Nelemans||2009| [Regimbau & Hughes||2010| [Regim-j 
jbau et al.|20T^ this is shown in Eig.j^with a green dashed line 

( jSathyaprakash & Schutz|2009^ . The duty-cycle for these events, 
for future space-borne gravitational observatories like DECIGO, 
will be higher than unity, i.e. the signals overlap and create a con¬ 
fusion noise of unresolved sources. See jYagi & Seto|p011| ( for a 
detailed description of the binary confusion noise as observed by 
the future DECIGO observatory. 

The gravitational signal generated by DCBHs is characterized 
by D <C 1, so there is a superposition of a discontinuous signal 
over a background noise: in this case, the component separation is 
easier to perform. Our Eig.|^is a simple proof of concept of a sig¬ 
nal processing method that allows to: (i) probe the presence of an 
underlying signal buried into a background noise and (ii) estimate 
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Figure 7. Simple example of a signal processing method to separate a signal 
with <C 1 from the confusion noise. The sum h(t) = hs(t) + hn{t), 
where hs{t) is the shot-noise signal and h„{t) is the background noise 
component, is shown on the bottom, in red. The blue circles and the green 
crosses are the auto-correlation function for h{t) and for h„ (t) only, respec¬ 
tively. The duty-cycle of the signal is X> = 0.2 (this value has been chosen 
for an easier visualization on the plot), while the population of collapsing 
DCBHs has a realistic distribution of masses (10'^ Mq < M < 10® Mq, 
spin paramenters (0.01 < a < 0.5) and densities (10® gcm~® < p < 
10® g cm“®). See the text for a complete description. 


its periodicity. The red solid line is the sum of hs{t), the signal 
from DCBHs collapse with V — 0.2 (this value has been chosen 
for an easier visualization on the plot), and the binary confusion 
background hn{t), as in |Regimbau & Hughes| ( |2()10| > but with ob¬ 
servations extended up to 2 = 15. The signal is completely buried 
into the noise. Nonetheless, the blue points are the auto-correlation 
function of the signal, which shows a clear structure. It starts from 
a value Aq ~ 0.4 and drops approximately after ~ 4 cycles of the 
signal to a value ~ Ao/2Q\ this time separation is highlighted by 
the vertical orange line. In addition, the periodic behavior of the 
auto-correlation function resembles the intrinsic oscillations of the 
ringdown phase during the collapse. As a comparison, the green 
crosses represent the auto-correlation function for the noise com¬ 
ponent only: the absence of any structure is evident. The two pa¬ 
rameters which greatly affect the appearance of the auto-correlation 
function are the mass M and the mean density p. The total duration 
of the GWs emission depends on: (i) the duration of the collapse 
phase (proportional to 1/ ^fp) and (ii) the duration of the ringdown 
phase, which depends on M (see the expression for vq in Sec.[^. In 
addition, the mass greatly affects the overall normalization (magni¬ 
tude) of the signal, see Fig.[^ For masses M > 10® Mq, the signal 
would not be buried into the fore/background, then its detection 
would be straightforward. If we allow the mass of the population 
to vary below M = 10® Mq and the density to vary as well, what 
happens is that the auto-correlation function may have several rises 
and falls, around the typical time separations between these signals, 
as clearly shown in Fig.|^ 

This simple method proves the existence of an underlying pe¬ 
riodic signal and provides a range of possible values for the pe¬ 
riod, i.e. the time separation needed for the auto-correlation func¬ 
tion to drop. Anyway, a possible future detection of this signal with 
Ultimate-DECIGO would require the application of more sophisti¬ 
cated algorithms for data analysis, similar to those that have been 
proposed for the LISA experiment (e.g. [Crowder & Comish|2007^. 


4 CONCLUSIONS 


4|. We 
Thigh, 


In this work, we have employed modem waveforms and the im¬ 
proved knowledge on the DCBHs formation rate to estimate the 
gravitational signal emitted by the formation of DCBHs at 13 < 
2 < 20, in the so-called “DCBHs Era” (see Yue et al.|2014| 
have thoroughly discussed the reasons why the, unexpectedly 1 
DCBHs formation rate is not in tension with current estimates of 
the local BHs mass density. Indeed, there are currently no observa¬ 
tional and/or theoretical indications ruling it out and further work is 
then required in this field, as also recent works (e.g. [Kormendy &| 
|Ho|2013{|Comastri et al.|2015^ suggest. 

We have investigated the effects of a wide range of masses, 
spin parameters and mean mass densities on the gravitational sig¬ 
nal, in order to build up a realistic ensemble of these high-redshift 
sources. We have shown that, despite the very high rate of events 
for all-sky surveys (~ 10^ yr~^ up to 2 ~ 20), the integrated sig¬ 
nal from these sources is characterized by a very low duty-cycle 


(D 


10 '^), i.e. it is a shot-noise. This is a consequence of 


the small duration of the GW emission from this kind of sources 
(ATj™ ~2-30s). 

Our results show that the gravitational signal lies above the 
foreseen sensitivity of the Ultimate-DECIGO observatory in the 
frequency range (0.8 — 300) mHz, with a peak amplitude Ogm = 
1.1 X 10“®“^ at Vmax = 0.9 mHz and a peak Signal-to-Noise Ratio 
of ~ 22 at!/ = 20 mHz. This amplitude is lower than the Galactic 
confusion noise in the same frequency band, generated by binary 
systems of compact objects. 

For a gravitational signal characterized by D <C 1 and buried 
into this confusion noise, we have provided a very simple signal 
processing method to prove the existence of underlying periodic 
oscillations such as those we expect from DCBHs. However, more 
advanced techniques will be required to separate this signal from 
background and foreground noise components. The signal inves¬ 
tigated in this paper lies in the same frequency range of the one 
produced by NS binaries and, in the same way, acts as a noise for 
the cosmological background, which is the ultimate target of the 
experiment DECIGO. For this reason, it is very important to model 
any gravitational signal that may fall into the frequency range of 
interest for the primordial signal. 

Despite all the technical difficulties, the actual observation of 
DCBHs may be a keystone in modem Cosmology, providing a 
significant contribution to the formation theory of Super-Massive 
Black Holes and to the understanding of the Cosmic Infrared Back¬ 
ground. 
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